function f=mynorminv(p,mu,sigma)

z=-sqrt(2)*erfcinv(2*p);
f=z*sigma + mu;

    %
